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RESEARCH MEMORANDUM 


EXPERIMENTAL DETERMINATION OF LINEAR DYNAMICS OF 
TWO -SPOOL TURBOJET ENGINES 
By David Novik and Herbert Heppler 


SUMMARY 

An experimental determination of the dynamics of two-spool turbojet 
engines is presented. Step inputs in fuel flow and in exhaust-nozzle 
area were utilized for initiation of transients on a test engine) and 
the resultant transient data were analyzed with fuel flow, turbine-inlet 
temperature, and exhaust-nozzle area considered as independent variables. 
Transfer functions, descriptive of general response, were obtained for the 
more important engine variables . 

From the test-program data, the speed responses of both spools to 
changes in turbine-inlet temperature appeared to be first-order lag 
systems with the same time constants. Data presented for another two- 
spool turbojet engine of different manufacture indicated that the speed 
responses of the two spools to changes in fuel flow also were identical 
first-order lags* 

With respect to exhaust -nozzle area changes, it was found that the 
outer-spool response appeared to be a first-order lag. Response of the 
inner spool near rated speed was apparently second order, but the magni- 
tude of change in the operating region was so small as to have no appre- 
ciable effect on the first-order outer-spool response or on the response 
of other engine variables . 


Comparison of engine response to the independent variables, turbine- 
inlet temperature (or fuel flow) and exhaust-nozzle area, indicated that 
engine dynamics for changes in turbine-inlet temperature (or fuel flow) 
were functions of inner-spool response and that engine dynamics for 
changes in exhaust-nozzle area were functions of outer-spool response* 

The engine responded more rapidly to a change in exhaust-nozzle area 
than to a change in turbine-inlet temperature or fuel flow* Experimental 
speed responses were found to agree with analyticallyr'derived equations 
of response . 
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INTRODUCTION 

In general, determination of the dynamics of an engine, supplemented 
by knowledge of the steady-state engine performance, is considered to he 
a prerequisite to the successful design of an optimum engine control 
system. Linear dynamics are considered necessary for stability and 
control-parameter considerations, and nonlinear dynamics are vital for 
delineation of surge limits and maximum acceleration potential. 

Sufficient effort has already been expended on the determination 
of single-spool turbojet- and turboprop -engine dynamics so that these 
engines are generally accepted as first-order lag systems whose responses 
can be characterized by means of time constants. The addition of the 
two-spool engine to those powerplants currently available for military 
aircraft use introduces an engine that contains an energy storage element 
with each spool and hence is inherently at least a second-order system. 
Published information pertinent to the dynamics of the two-spool engine 
has thus far been relatively scarce. 

At the NACA Lewis laboratory, experimental programs have been under- 
taken for the purpose of obtaining both linear and nonlinear dynamics of 
two-spool engines, in conjunction with steady-state performance under 
varying conditions of operation. The scope of this report is confined 
to an investigation of the linear dynamics of two-spool engines with the 
expectation that the information presented will be useful in control 
synthesis from stability and control-parameter considerations . A two- 
spool engine, operated in a sea-level test stand, was subjected to step 
inputs in fuel flow and exhaust-nozzle area. Some of the original tran- 
sient data obtained are presented, together with typical results of 
harmonic analysis of the data and the transfer functions of the engine 
variables derived therefrom. An acceleration trace for a two- spool 
engine of different manufacture is also shown and analyzed. 


SYMBOLS 

The following symbols are used in this report: 
A exhaust-nozzle area 

a rise ratio 

F thrust 

I moment of inertia 

K constant 
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f 

y f 

y t 


N engine speed 

P total pressure 

p differential operator 

T total temperature 

t time 

W.p fuel flow 

engine variable at end of transient 
engine variable at time equals t 
y-Q engine variable at time equals zero 

6 pressure correction, ratio of ambient pressure to NACA standard 

sea-level pressure 

0 temperature correction, ratio of ambient temperature to NACA 

standard sea-level temperature 

X time constant with T^ and A considered as independent 

variables 

x ' time constant with and A considered as independent 

variables 

cp phase shift, deg 

os frequency, radians/sec 

Subscripts : 

1 inner 

0 outer 

1 outer-compressor inlet 

2 outer-compressor discharge and inner-compressor inlet 

3 inner-compressor- discharge 

4 turbine inlet 

9 tail pipe 
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Functions, transfer functions, and partial derivatives: 

Q inner-spool torque expressed as function of N , Nj_, T^, and A 

inner-spool torque expressed as function of N Q , N i , W f , and A 

outer-spool torque expressed as function of N Q , N^_, T^, and A 

outer-spool torque expressed as function of N Q , N-j_, W^, and A 


Q' 

<1 


Aq_ 

AT, 


transfer function for response of outer-spool torque to a change 
in turhine-inlet temperature at constant exhaust-nozzle area 


% 




* (st), 

V ol 


i 


APPARATUS t 
Engine 

The two-spool engine used in this investigation is presented sche- 
matically in figure 1. The engine is essentially composed of an outer 
axial-flow compressor connected hy a through shaft to a turbine , an 
inner compressor connected to a turbine,, and a multiunit semiannular 
combustion chamber. Automatic bleed valves, closed during high-speed 
operation and open at low speeds, are located between the two com- 
pressors. Experimentally determined steady-state performance curves 
are shown in figure 2 as functions of inner-spool speed. 


Fuel System 

The engine fuel system, up to the flow divider, was supplanted by 
a test-facility fuel system as shown in figure 3. The fuel control con- 
sisted of a differential relief valve controlling the pressure drop 
across a throttle valve, the position of which (and thus the area) was 
varied by an electrohydraulic servomotor. In such a system, fuel flow 
is proportional to throttle -valve position. The natural frequency of 
the differential relief valve was approximately 250 cycles per second 
(ref. l), and the natural frequency of the electrohydraulic servomotor 
was approximately 60 cycles per second. Because of the high natural 
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frequency of the differential relief valve, the fuel flow was propor- 
tional to the throttle -valve position for all throttle -valve transients. 
Therefore, throttle-valve position was used as the recorded indication 
of fuel flow* The over-all response of the fuel system, up to the engine 
flow divider, was underdamped hut reached equilibrium within 0*03 second 
following a step input* 


Variable Exhaust-Nozzle Area 

Exhaust-nozzle area was varied by means of pneumatically operated 
flaps as shown in figure 4, The two flaps were linked together so that 
they moved simultaneously, and motion of the flaps from open to closed 
approximated a ramp of about 0*04 to 0*05 second* 

Engine speeds * - Transients in engine speeds were measured by re- 
cording the d-c voltage outputs of electronic counters that counted the 
pulses supplied from 180-tooth tachometer generators . Before recording, 
the speed signals were filtered with resistance.- inductance-capacitance 
filters that had characteristics as shown in figure 5* 

Engine temperatures * - Temperatures at the discharge of the outer 
and inner compressors were measured by three iron- constant an thei uncouples 
in series at each station* Thermocouples were made of 27-gage wire and 
ground to a sharp arrowhead point for fast response without compensation* 

Turbine-inlet and tail-pipe temperatures were measured by three 
chromel-alumel thermocouples in series at each station. These high- 
temperature thermocouples were made of 14-gage wire and were compensated* 
The compensated signal to the recorder had a lag of approximately 0.03 
second. 

Engine pressures . - Engine pressures were measured by means of 
strain-gage pressure pickups. Pressure responses had apparent time 
constants below 0.01 second. 

Thrust . - Engine thrust was measured by a strain-gage thrust link. 


EXPERIMENTAL PROCEDURE AND DATA 
Test Procedure 

Data indicative of the response of outer-spool speed, inner-spool 
speed, outer-compressor-discharge pressure. and temperature, inner- 
compressor-discharge pressure and temperature, turbine-inlet tempera- 
ture, tail-pipe pressure, tail-pipe temperature, and thrust were re- 
corded during transients induced by changes in fuel flow and exhaust- 
nozzle area. 
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Fuel-flow transients . - Fuel-flow transients, covering a range of 
engine speeds, were initiated from approximate step inputs in fuel flow. 
Fuel changes were of such a magnitude as to result in inner- spool speed 
changes of approximately 3 percent on the assumption that these speed 
changes were small enough to insure linearity. Exhaust -nozzle area was 
fixed during fuel-flow transients . . 

Exhaust-nozzle-area transients . - Exhaust- nozzle -area transients 
were induced hy closing the flaps installed on the exhaust nozzle. 
Closing the flaps corresponded to a change in exhaust-nozzle area from 
105 percent area to 100 percent area and resulted in outer-spool speed 
changes of about 3 percent. Fuel flow was maintained constant during 
exhaust-nozzle-area transients. 


Transient Data 

Fuel-flow transient data . - Samples of transient data obtained for 
approximate step inputs in fuel flow are shown ih figures 6(a) and (b) . 

It is of interest to note that both spools appear to move simultaneously 
as though the rotors were geared together. An increase in response time 
is evident for the transient at the lower engine speed. 

Figure 6(c) has been included to illustrate the discontinuities in 
the engine variables when the compressor bleed valves automatically close 
during a transient . 

Exhaust-nozzle-area transient data . - Typical transient data ob- 
tained for approximate steps (ramps of 0.04 to 0.05 sec) in exhaust- 
nozzle area are shown in figure 7. In. all the transients shown, the 
exhaust nozzle went from an open to a closed position. Area traces 
may be erratic or are not recorded at all because of a faulty potenti- 
ometer coupling. 

Inspection of exhaust-nozzle-area transients reveals that the inner- 
spool response is of extremely small magnitude except at the low engine 
speeds and that its net change near top speeds is opposite in direction 
from that of the outer spool. Net changes in thrust for changes in 
exhaust-nozzle area are always small, despite the initial response. 


ANALYSIS OF DATA 
Theoretical Analysis 

In order to facilitate evaluation of the experimental data obtained, 
reference is made to the following equations of response. These equations 
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are derived in appendix A for engine operation near rated speed. 
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Although equations (A13) to (A16) specify turhine-inlet temperature 
rather than fuel flow as an independent variable, these equations are 
also indicative of fuel-flow response near rated speed because turbine- 
inlet temperature response to fuel flow normally approaches unity near 
rated engine speed. 


Response of Engine Variables to Fuel Flow from Harmonic 
Analysis of Transient Data 

Typical frequency-response curves obtained from harmonic or Fourier 
analysis (ref* 2) of fuel step transients are shown in figures 8 and 9. 
Frequencies corresponding to the points shown are given in table I* 
Cursory examination of these frequency-response curves shows higher- 
order responses that are obviously not in agreement with the responses 
anticipated from analysis. For example, it can be seen that the inner- 
spool speed response, even at high speed,, appears to be considerably 
more complex than the anticipated approximation of first-order response* 
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It is "believed that the apparent discrepancy "between analysis and 
data can he attributed to incomplete combustion during initial phases 
of the transients, and possibly to dynamics of the spring-loaded flow 
divider. Combustion dynamics are implied from inspection of figure 8(g), 
which shows the response of turbine-inlet temperature to fuel flow at 
96.1 percent rated inner-spool speed. The large phase shifts evident at 
high frequencies cannot be attributed to engine-speed dynamics because 
the speed responses are very much attenuated at the high frequencies . 

The higher-order phase shifts in the temperature response must therefore 
be indicative of combustion lags immediately following the fuel step. 

Despite the combustion effects, some general trends in the responses 
to fuel flow can be ascertained. The frequency-response curves of the 
two speeds (figs. 8(a) and (b) or 9(a) and (b)) are essentially the same, 
indicating that one transfer function defines the response of both spools. 
Smaller phase shifts for other engine variables are indicative of leads 
in the response of these variables as compared with the speed response. 
Engine response is slower at the lower engine speeds as evidenced by the 
increased phase shift for a given frequency at the lower engine speed. 

The transient data of figure 10, taken from a two-spool turbojet 
engine of different manufacture than the test engine, can be analyzed 
to give more valid indications of response to fuel flow because of the 
slower rate of fuel addition. The frequency-response curves of outer- 
and inner-spool speeds with respect to fuel flow for these data are 
shown in figure 11. It can be seen that the responses of both spools 
appear to be identical first-order lags. The response of the inner 
spool is therefore in agreement with the response anticipated from 
equation (A14) . The response of the outer spool, however, is simpler 
than anticipated by equation (A13) in that it is the same first-order 
lag as that of the inner spool. 


Response of, Engine Variables to Turbine-Inlet Temperature 
from Harmonic Analysis of Transient Data 

Inasmuch as the turbine-inlet temperature reflects the amount of 
energy actually available for acceleration, the frequency response of 
engine variables with respect to turbine-inlet temperature can be ex- 
pected to yield more predictable results than the responses to fuel- 
flow changes . 

Typical frequency-response curves showing the responses of engine 
variables to turbine-inlet temperature (obtained from harmonic analysis 
of the fuel-step transient data with reference to turbine-inlet tempera- 
ture as the input) are shown in figures 12 and 13. 
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Speed responses . - The frequency responses of outer-spool speed to 
turbine- inlet temperature are shown in figures 12(a) and 13(a ) } and the 
responses of inner-spool speed are shown in figures 12(h) and 13(h) . 

The responses of hoth spools appear to he first order # and the indicated 

time constants {% = i ^ of each spool are 

frequency at 45° phase shift 

almost identical. The response of hoth spools to turhine-inlet tempera- 
ture appears to he 


AN _ K 
AT 4 Tp + 1 


where K represents the steady-state gain of the spool being considered 
and t represents a time constant common to hoth spools .. 


The experimental response of inner-spool speed to changes in turhine- 
inlet temperature therefore agrees with the analysis (eq. (A14)) in that 

“ I i 

it is a first-order response with a time constant T-? = — — . The fact 


that the experimental response of outer-spool speed appears to he virtually 
the same as the response of inner-spool speed means that the second-order 
system of equation (A13) must degenerate into a first-order system with 
a time constant It is assumed that this simplification occurs in 

the following manner: 


AN q 

AT. 



(t ± p + l) + 
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% (^iP + D(t o P + l) 


(A13) 


■but apparently 

(t .p + l) + cx t q P + 1 

1 ^^4 


and therefore 

K| .. % i 


Inasmuch as figure 11 has already indicated similar responses for the two 
spools of an entirely different engine^ the preceding factoring may he 
considered more of a generality than a coincidence. 
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A plot of the variation of the inner-spool time constant as a 

function of inner- spool speed (measured at the end of each transient) 

is shorn, in figure 14.* Variation of the apparent time constants T ob- 
tained for the outer-spool response to turbine -inlet temperature is also 
shown* Figure 14 shows no appreciable difference between the two re- 
sponses although there appears to be a very slight tendency for the 
inner spool to be slower than the outer spool at high speeds* Compressor 
bleed at 92*5 percent rated inner-spool speed is seen to reduce the time 
constant by an almost negligible amount* The over-all trend of both 
spools exhibits the usual variation of response time with engine speed 
normally associated with single-spool turbojet engines* 


Outer-compressor-discharge pressure response * - Response of outer- 
compressor-discharge pressure P£> to changes in turbine-inlet tempera- 
ture at constant area (figs* 12(c) and 13(c)) appears to have a small 
lead with respect to the speed response. In transfer function form, 


K(arr iP + l) + X ) 


T iP + 1 


T iP + X 


* * Po 

The value of a (or — _) is the rise ratio obtained from the frequency- 


response curve such that a = - — — = an T ) li'fcu-de Infinite frequency # 

i amplitude at zero frequency 

(Where higher-order lags appear, they are neglected and the amplitude 
at infinite frequency is taken at the point of intersection of the 
faired curve with the zero phase-shift axis.) Variation of rise ratio 
with engine speed is shown in figure 15, from which it may be seen that 
the rise ratios are small and disappear at low engine speeds. 


Outer-compressor-discharge temperature response * - The outer- 
compressor-discharge temperature T^ exhibits a small degree of lead 
as shown in figures 12(d) and 13(d)* Data were too inconclusive to 
permit a plot of rise ratios) however, a rise ratio of approximately 
0.1 is estimated at rated speed. It is assumed that the trend of rise 
ratios with engine speed is similar to that shown for outer-compressor- 
discharge pressure (fig. 15)* 


Inner-compressor-discharge pressure response . - Figures 12(e) and 
13(e) indicate that the response of inner-compressor-discharge pressure 
P3 to turbine-inlet temperature at constant area is a lead-lag system 
that can be approximated by the relation 

AP 3 _ K(aT iP + 1) _ K fv + X ) 

AT4 ^ T iP + 1 ^ TT iP + 1 
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Values of rise ratio for inner-compressor-discharge pressure as a func- 
tion of engine speed are shown in figure 16, from which it can he seen 
that the rise ratio is small and increases with increasing engine speed. 

Inner-compressor-discharge temperature response , - The response of 
inner-compressor-discharge temperature to turbine-inlet temperature 
appears to he a lead-lag system (figs. 12(f) and 15(f)) as would he ex- 
pected from the lead-lag response obtained for inner-compressor-discharge 
pressure. The’ temperature rise ratios are plotted in figure 17 and ap- 
parently are quite small with a very slight increase at increasing engine 
speed c 

Tail-pipe pressure response . - Figures 12(g) and 13(g) show that 
the tail-pipe pressure response to a change in turbine-inlet tempera- 
ture appears to he a lead-lag system with small rise ratios. As shown 
in figure 18, the rise ratio increases with increasing engine speed* 

Tail-pipe temperature response , - Rise ratios indicative of tail- 
pipe temperature lead-lag response can he obtained from figures 12(h) 
and 13(h), If the premise is accepted that combustion within the com- 
bustors was incomplete, then roughness of the data in the frequency- 
response curves is a reasonable result. Large rise ratios that de- 
creased with Increasing engine speed were obtained as shown in figure 19. 

Thrust response , - Thrust response to turbine-inlet temperature 
( figs • 12 ( i ) and 13 ( i ) ) appears to be lead-lag. Variation of rise ratio 
with engine speed is shown in figure 20* Rise ratios are small and 
increase with increasing engine speed. 

Effect of large fuel steps on appearance of lead response , - In 
order to approach linearity, small steps in fuel flow were used to Initiate 
the transients. The over-all changes in engine variables were therefore 
small, and the leads or immediate changes inherent in the dynamics of 
some engine variables are not visibly apparent in the transient data. 

The leads become more obvious for large fuel steps, as shown in figure 21, 
in which the leads in turbine-inlet temperature, tail-pipe temperature, 
thrust, and inner-compressor-discharge pressure can be clearly seen. 


Response of Engine Variables to Changes in Exhaust-Nozzle Area from 
Harmonic Analysis of Transient Data and from Semilog Plots 

Illustrative frequency-response curves showing the response of engine 
variables to steps in exhaust-nozzle area at constant fuel flow are shown 
In figures 22 to 25, These curves were obtained from harmonic analysis 
of the transient data shown in figure 7 . Inspection of these frequency- 
response curves indicates characteristic responses that appear to be very 
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close to first order. Inasmuch as first-order exponential curves plot 
as straight lines on semilog paper , the transients of figure 7 have been 
normalized according to figure 26 and plotted on semilog paper in figures 
27 to 30, Scatter of. data toward the end of the semilog plots' is con- 
sidered to he the result of engine drift and difficulty of measurement, 
in addition to small higher-order engine dynamics . 

Outer-spool speed response . - Response of outer-spool speed to 
changes in exhaust -nozzle area at constant fuel flow is characterized 
by a 90° phase-shift semicircle in the frequency-response plots (figs. 
22(a), 23(a), 24(a), and 25(a)) and by a straight line on the semilog 
plots (figs. 27 to 30). Both methods of data analysis are therefore 
indicative of first-order response and tend to corroborate the response 
anticipated from equation (A15) . This response is now at constant' fuel 
flow rather than constant turbine-inlet temperature: 


where x = 

o q * 

H N 

o 

A plot of time constants for the response of outer-spool speed to 
exhaust-nozzle area, as obtained from frequency- response curves and 
from semilog plots, is shown in figure 31, as a function of engine speed. 
The time constants decrease as engine speed increases, and the response 
is essentially constant in the top-speed range. 


AN„ 


AA 


K 


W-f 


oc 


T^p + 1 


Inner-spool speed response . - Inner-spool speed response to changes 
in exhaust-nozzle area is considerably different from the response of 
outer-spool speed. As noted previously in the discussion of the transient 
data, the steady-state gain or net change of the inner spool is in the 
same direction as that of the outer spool at .low speeds (speeds below 
normal operating range), is rather indecisive and small at about 85 to 
90 percent of rated outer-spool speed (fig. 7(b)), and is small and op- 
posite in direction at speeds approaching rated speed (fig. 7(a)). 


The semilog plot (fig. 27) and the frequency-response curve (fig. 
22(b)) for the transient most closely approaching rated speed indicate 
that the inner-spool speed response is second order near rated speed. 

It is surmised that the transfer function of inner-spool speed is (from 


eq. (A16)) 


ANi 
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W* 
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(o-.’p + l) (t 'p + 1) 
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Although the inner spool has a dynamic response to changes in exhaust 
nozzle area^ as approximated in the preceding equation* the steady-state 
gain or net change is so small as to leave the outer-spool speed essen- 
tially unaffected and* to all outward appearances* first order. 

From near rated speed to the vicinity of 65 percent rated speed 
(based on outer-spool speed) * there is a gap in which the change in 
inner-spool speed with exhaust-nozzle area is small and inner-spool 
dynamics are essentially negligible and indeterminate. This region is 
typified by the transient of figure 7(b) and by the responses shown in 
figures 23 and 28. The essentially first-order appearance of outer- 
spool speed indicates that the change in inner-spool speed is too small 
to have any measurable effect. 

In the vicinity of 65 percent rated outer-spool speed* the inner- 
spool response to exhaust -nozzle area appears first order with a longer 
time constant than that of the outer spool (figs. 24(b) and 29). The 
response is large enough to have a small effect on the outer-spool speed* 
toward the end of the transient* as indicated by' the break in slope of 
the outer-spool speed in figure 29. 

At the lowest speed for which transient data were obtained (fis. 
7(d))* corresponding to a final outer-spool speed of 44.7 percent of 
rated* the inner-spool response continued to appear first order (figs. 
25(b) and 30) . The magnitude of change has increased to approximately 
the same as that of the outer spool* and the time constant has become 
identical to that of the outer spool . 

Basic engine response to change in exhaust-nozzle area . - Plots of 
engine variables (other than speed) in figures 27 to 30 fall parallel to 
the outer-spool speed traces* and in figures 22 to 24 the maximum ampli- 
tudes of the engine variables occur at frequencies that correspond to 
the outer-spool time constants. The inner-spool dynamics therefore have 
essentially no effect on the dynamics of the engine variables* and rise 
ratios can be obtained on the premise that the first-order lag of the 
outer spool represents the basic engine response. 

Outer-compressor-discharge pressure response . - The frequency- 
response curves of figures 22(c)* 23(b)* and 24(c) show no definite 
dynamics in outer-corapressor-discharge pressure response other than the 
first-order speed lag. The semilog plots of figures 27 to 30 show in- 
dications of a slight lead with respect to speed* but it is believed that 
this is because compensation was not made for the lags in the speed traces 
It is concluded that the response of outer-compressor-discharge pressure 
to changes in exhaust-nozzle area is basically a first-order lag that 
follows the change in outer-spool speed. 
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Outer-compressor-discharge temperature response . - Response of 
outer-compressor-discharge temperature to changes in exhaust-nozzle area 
may he seen from the frequency-response plots of figures 22(d) and 23(c) 
and from the semilog curves of figures 27 and 28. As in the response 
of the outer-compressor pressure, the temperature appeared to follow the 
outer-spool speed without additional dynamics . Transient data for outer- 
corapressor-discharge temperature were subject to large fluctuations such 
that analysis in some cases was impossible, and the faired curves used 
in analysis of the data must be construed as rough approximations . 

Inner-compressor-discharge pressure response . - Response of inner- 
compressor-discharge pressure to changes in exhaust-nozzle area may be 
seen from the frequency-response curves of figures 22(e), 23(d), and 
24(d) and from the semilog plots of figures 27 to 30. The response 
apparently follows that of the outer spool without any other dynamics, 
despite inner-spool speed dynamics. 

Inner-compressor-discharge temperature response . - The original 
transient-data traces of inner-compressor-discharge temperature were 
subject to large fluctuations such that some traces could not be faired 
at all for analysis and other faired traces were rough approximations . 
Typical analyzed responses obtained are shown in the frequency-response 
curves of figures 22(f) and in the semilog plot of figure 27. The 
response of inner-compressor-discharge temperature appears to be the 
same as that of the outer-spool speed. 

Turbine-inlet temperature response . - The frequency-response curves 
of turbine-inlet temperature (figs. 22(g), 23(e), and 24(e)) and the 
semilog plots of figures 27 to 29 indicate that turbine-inlet tempera- 
ture response to exhaust-nozzle area is a lead-lag system. 

The rise ratios (indicative of initial jump) are very small and 
approach zero in the top-speed range. The small values of rise ratio 
coupled with the fact that the steady-state change with area is small 
implies that for small changes in exhaust-nozzle area the turbine-inlet 
temperature response could probably be assumed to be zero. Rise ratios 
obtained are plotted in figure 32 . 

Tail-pipe pressure response . - Response of tail-pipe pressure to 
exhaust-nozzle area may be seen from the frequency-response curves of 
figures 22(h), 23(f), and 24(f) and from the semilog plots of figures 
27 to 30. Tail-pipe pressure response has a very definite lead as ex- 
pected for exhaust-nozzle area changes. The second-order lags notice- 
able at the high-frequency end of some of the frequency-response curves 
may possibly be attributed to inner-spool speed changes. If these second- 
order effects are neglected, a plot of rise ratios for tail-pipe pres- 
sure as a function of engine speed is obtained as shown in figure 33. No 
discernible trend in rise ratio was found. 
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Tail-pipe temperature response , - Tail-pipe temperature response to 
changes in exhaust -nozzle area is indicated hy the frequency-response 
curves of figures 22(i), 23(g ) , and 24(g) and hy the semilog plots of 
figures 27 to 29, A small amount of lead is apparent in the response 
and the lag decay follows outer-spool speed, 

A first-order lead-lag system is used as an approximation of the 
response, and rise ratios are small and increase with increasing engine 
speed. The estimated point shown in figure 34 is the result of allowance 
made for the fact that electrical compensation was not used to correct 
for thermocouple lag in the tail-pipe temperature measurement of the 
transient shown in figure 7(a). Inasmuch as the thermocouple time con- 
stant could he reasonably estimated at 0,3 second and yet the phase 
shift resulting from this 0.3-second lag was more than cancelled out hy 
the lead of the calculated response (fig. 22(i)); it was assumed that 
the lead time constant must have heen at least 0,3 second. 

Thrust response . - Frequency-response curves for the response of 
thrust to changes in exhaust -nozzle area are shown in figures 22(j), 

23(h ) } and 24(h ) , and relative slopes are shown in the semilog plots of 
figures 27 to 29. If positive rise ratios are assumed for the condition 
in which the final value of thrust exceeds the initial value of thrust 
for a reduction in exhaust -nozzle area, then the rise ratios were found 
to he positive at the low engine speeds and negative at the high engine 
speeds, A plot of rise ratio as a function of engine speed is shown in 
figure 35 , Infinite values of rise ratio indicate a condition in which 
initial and final values of thrust are identical. The initial jumps in 
thrust^ corresponding to 5~percent exhaust-nozzle area changes , were on 
the order of 5 percent. 


Transfer Functions 

Responses of all engine variables Investigated have heen approxi- 
mated as transfer functions and are listed in table II. The transfer 
functions are based on the inner-spool time constant for fuel flow and 
turbine-inlet temperature inputs and are based on the outer-spool time 
constant for exhaust-nozzle area inputs . 


Relative Response 

Comparison of the time-constant curves of figures 14 and 32 shows 
that the engine responds more rapidly to a change in exhaust-nozzle 
area than it does to a change in turbine-inlet temperature , Wear rated 
speed,, the response of turbine-inlet temperature to fuel flow approaches 
unity (fig. 8(g)) and; hence ; engine response to exhaust-nozzle area is 
also faster than the response to fuel flow, at least near rated speed. 
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It may also "be concluded that, if each spool is taken as a separate 
entity, the outer spool has a smaller time constant (or energy storage 
element) than the inner spool. This is deduced from the fact that engine 
response to changes in exhaust -nozzle area is a function of outer-spool 
time constant, whereas engine response to changes in fuel flow appears 
to be a function of inner-spool time constant . 


SUMMARY OF RESULTS 

Transient data obtained from step inputs in fuel flow and exhaust - 
nozzle area for sea-level operation of a two-spool turbojet engine have 
been presented and analyzed. Because the fuel steps resulted in a non- 
linear time variable relation between fuel flow and turbine-inlet tempera- 
ture, the effect of fuel-flow transients on engine dynamics was analyzed 
with respect to turbine-inlet temperature as the input. The speed re- 
sponses of each spool to changes in turbine-inlet temperature were es- 
sentially identical and could be characterized as a first-order lag. 

Time constants decreased with increasing engine speed in a manner similar 
to the time-constant variation of a single-spool turbojet engine. Data 
obtained from an entirely different two-spool engine were also presented 
in order to determine the response to fuel flow and to permit greater 
generalization of results . It was found that the responses of both spools 
to a change in 'fuel flow again appeared to be identical first-order lags. 

The response of engine variables to a change in exhaust -nozzle area 
at constant fuel flow was found to be more complicated than the response 
to turbine-inlet temperature or fuel flow. It was found that the outer- 
spool response appeared to approximate a first-order lag, but that the 
inner-spool response appeared to be higher order in the operating speed 
range. However, the magnitude of the inner-spool response was extremely 
small and consequently did not appreciably affect the apparent first-order 
response of the outer spool. At low speeds the magnitude of inner-spool 
speed changes increased and the response became first order. 

As a general rule it appeared that engine dynamics were functions of 
inner-spool response for changes in turbine-inlet temperature (or fuel 
flow) and were functions of outer-spool response for changes in exhaust- 
nozzle area. Response to exhaust -nozzle area was faster than response 
to fuel flow or turbine-inlet temperature. 

Experimentally determined responses of inner-spool speed to fuel 
flow and to turbine -inlet temperature and of both inner- and outer- 
spool speeds to exhaust -nozzle area were found to agree With analytically 
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derived equations of response. The conclusion that inner-spool time 
constant represents the basic engine response to fuel flow and turhine- 
inlet temperature was drawn from analysis. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland^ Ohio,, October 18,, 1954 
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APPENDIX - DERIVATION OF RESPONSE EQUATIONS 

The response of each spool to changes in the independent variables , 
turbine-inlet temperature and exhaust-nozzle area, is obtained from 
linearization of the functional relations: 


<1 = q(N 0 ,N t ,T 4 ,A) 
Q = Q(N 0 ,Ni,T 4 ,A) 


(Al) 

(A2) 


Linear expansion of the functional relations around steady-state operating 
points gives 

dN_ , , 

Aq = I Q A — I 0 pAN 0 = <i No AN 0 + q^AN.^ + q. T ^AT 4 + q^AA (A3) 

dN. 

AQ = I ± A — - = IiPANi = Qjj AN q + Q K A% + Qij AT 4 + Q^AA (A4) 

CLu O 1 ^ 

By dividing equation (A3) by equation (A4) and eliminating either 
AN^ or AN Q , equations for the response of each spool can be obtained. 

At constant area, AA = 0) and 


ANq 


q^^iP - + % 4 ( iN i 

at 4 

A 

( r iP - Q Ni )(I 0 P - qN 0 ) - ^Nc^Ni 

AN. 

l 


% (IqP " ) + % 

4 o o 4 

AT 4 

A 

(3-iP - % i )( I oP " < lN 0 ) “ 


(A5) 


(A6) 


Equations (A5) and (A6) can be rearranged as follows: 



B _ V V»i 

p + 1 j + - - 


p + 1 


O'No'lNi 

"o 


(A7) 
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AW i 


AT, 


Qm 




Q N 0 1T 4 

+ 


%-L Z" 1 ! \ /" -I o \ % 0 *Ni 

V% p + 1 Jv% p + 1 J' %% 


-I© -I-? 

Let t; = and t, = __ij then, 

° % %i 


(A8) 


AN. 


AT, 


q T, (^iP + l) + 014% 


ANi 


AT, 


A X 




(t^P + 1)(t o P + 1) - 


l N 
o XN i 


(^qP + l) + 


V*o 

%</*T 4 


|A \ V H i 

1 (TlP + I)(T ° P + x) * 5 a 


(A9) 


(A10) 


If equations (A3) and (A4) are used to determine response to exhaust 
nozzle area at constant turbine-inlet temperature (AT 4 = 0) , then the 
following equations are obtained in a manner similar to the preceding 
development : 

q a%^ 

(t,p + 1) + 


AN 

c 

AA 


%i q A 


T 4 q N 0 


(tiP + 1)(t g P + l) - 


% 0 ^N i 


(All) 


AN, 


AA 




Q N 0 q A 

(t ° p + p) + 


T„ 


\ q n o % 1 

Kp + 1)(T 0 P + 1) - 

i o 


(A12) 
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Evaluation of the partial derivatives on the basis of thermodynamic 
engine relations near rated-speed conditions indicates that the quantities 

— j are negligible, the quantity Q A is equal to zero, 


Q Ni% c 


and 




and the quantity 


^N^A 

V N o 


is very small. 


The equations of response for a two-spool engine near rated speed 
are therefore approximated by the following equations: 


I AT. 


(x.p + 1) + - 

_ % Q N i ( ^T 4 

A % (TjP + 1)(t o P + 1) 


(A13) 


an 4 


AT. 




A 


Qi 


'N 


■ T 1 p + 1 
i 1 


(A14) 


^o 

AA 


AW i 


AA 


= 1 
T 4 % q t oP + 1 

% 0 q A 

^N.^N 
1 o 


(t ± P + 1)(t q P + 1) 


0 


(A15) 


(A16) 
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TABLE I. - FREQUENCIES CORRESPONDING TO POINTS 
SHOWN ON FREQUENCY-RESPONSE CURVES 


Point 

Frequencies applicable for figs.— 


8, 11, 12, 21, 22 

9, 13, 23 

24 

A 

0.1047 

0.05236 

0.0262 

B 

.131 

.06545 

.0327 

C 

.157 

.07854 

.03295 

D 

.2094 

.1047 

.05325 

E 

.262 

.1309 

.0654 

F 

.314 

.1571 

.0785 

G 

.419 

.2094 

.1047 

H 

.523 

.2618 

.131 

I 

.628 

.3142 

.157 

J 

.873 

.4363 

.218 

K 

1.047 

.5263 

.262 

L 

1.31 

.6545 

.327 

M 

1.57 

.7854 

.3295 

N 

2.094 

1.047 

.5235 

0 

2.62 

1.309 

.654 

P 

3.14 

1.571 

.785 

Q 

4.19 

2.094 

1.047 

R 

5.23 

2.618 

1.31 

S 

6.28 

3.142 

1.57 

T 

8.73 

4.363 

2.18 

U 

10.47 

5.236 

2.62 

V 

13.11 

6.545 

3.27 

W 

15.71 

7.854 

3.925 

X 

20.94 

10.47 

5.235 
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Corrected Inner-spool speed, percent rated bleeds closed 


(a) Outer-spool speed, fuel flow, and thrust. 

Figure 2. - Steady-state performance of two-spool turbojet engine. 
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Corrected inner-spool speed, percent rated bleeds closed 


( b ) Temperatures . 


Figure 2. - Continued « Steady-state performance of two-spool turbojet 
engine , 
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Figure 2. - Concluded, 
engine . 
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Figure 3. - Fuel system. 
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Figure 5. - Characteristic filters used on engine-speed measurements. 
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(a) Corrected outer-spool speed range, 95.9 to 92.9 percent rated; corrected inner-spool speed range, 
96.6 to 96.7 percent rated; bleed valves closed. 



(b) Corrected outer-spool speed range, 88.0 to 85.1 percent rated; corrected inner-spool speed range, 
92.9 to 93*0 percent rated; bleed valves closed. 


Figure 7. - Response of engine variables to change in exhaust-nozzle area. 
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(c) Corrected outer-spool speed range/ 65*9 to 63.1 percent rated; corrected inner- spool speed range, 
83.8 to 83.1 percent rated; "bleed valves open. 



(d) Corrected outer-spool speed range, 46.4 to 44.7 percent rated; corrected inner-spool speed range, 
71.5 to 70.0 percent rated; "bleed valves open. 


Figure 7. - Concluded. Response of engine variables to change in exhaust-nozzle area. 
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(b) Inner-spool speed. 
Figure 11. - Frequency-response curv-es for engine of different manufacturer. 
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(f) Inner-compressor-discharge temperature. 



Phase angle, deg 
(g) Tail-pipe pressure. 



(h) Tail -pipe temperature. 


(i) Thrust. 


Figure 13. - Concluded. Frequency response of engine variables to change in turbine-inlet temperature for increase in corrected 
outer-spool speed from 57.6 to 62.9 percent rated and increase in corrected inner-spool speed from 79.6 to 82.5 percent rated. 
Bleed valves open. 
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rrected inner-spool speed, percent rated 


Figure 17. - Variation of inner- compressor- discharge- temperature rise ratio 
with engine speed for change in turbine- inlet temperature. 



Corrected inner- spool speedy percent rated 


Figure 18. - Variation of tail-pipe-pressure rise ratio with engine speed 
for change in turbine-inlet temperature. 
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Corrected inner-spool speed, percent rated 


Figure 19. - Variation of tail-pipe-temperature rise ratio with engine speed 
for change in turbine-inlet temperature. 
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Phase angle, deg 
fa) Outer-spool speed. 



-100 -90 -00 -70 -60 -50 -40 

Phase angle, deg 

(b) Inner-spool speed. 



( <3 ) Outer-compressor-diseharge pressure. (d) Outer-compressor-discharge temperature. 



Phase angle, deg 

(e) Inner-compressor-discharge pressure. 

Figure 22. - Frequency response of engine variables to change in exhaust -nozzle area 
for decrease in corrected outer- spool speed from 96.9 to 92.9 percent rated and 
increase in corrected inner- spool speed from 96.6 to 96-7 percent rated. Bleed 
valves closed. 
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(a) Outer-spool speed. (b) Inner-spool speed. 



(c) Outer-compressor-discharge pressure. (d) Inner-compressor-discharge pressure. 

Figure 24. - Frequency response of engine variables to change In exhaust-nozzle area for decrease in corrected outer-spool speed from 
65.9 to 63.1 percent rated and decrease in corrected inner-spool speed from 83.8 to 83.1 percent rated. Bleed valves open. 
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(a) Outer-spool speed. 

Figure 25. - Frequency response of engine variables to change in e: 
46.4 to 44.7 percent rated and decrease, in corrected inner-spool 
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Phase angle, deg 

(b) Inner-spool speed. 


-nozzle area for decrease in corrected outer-spool speed from 
from 71.5 to 70.0 percent rated. Bleed valves open. 
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Outer- spool speed 
Inner- spool speed 
Turbine- inlet temperature 
Tall- pipe temperature 
Outer- compressor- discharge 
pressure 

Inner- compressor- discharge 
pressure 

Tail-pipe pressure 
Outer- compressor- discharge 
temperature 

Inner-compressor- discharge 
temperature 
Thrust 


Figure 27. - Semilog plot of response of engine variables to change in 
exhaust- nozzle area for decrease in corrected outer- spool speed from 
95.9 to 92.9 percent rated and increase in corrected inner- spool speed 
from 96.6 to 96.7 percent rated. Bleed valves closed. 
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Outer- spool speed 
Turbine- inlet temperature 
Tail- pipe temperature 
Outer- compressor- dis charge 
pressure 

Inner- compressor- dis charge 
pressure 

Tail- pipe pressure 
Outer-compressor-discharge 
temperature 
Thrust 


3 12 3 

Time, sec 

Figure 28. - Semilog plot of response of engine variables to 
change in exhaust -nozzle area for decrease in corrected 
outer- spool speed from 88.0 to 85.1 percent rated and in- 
crease in corrected inner- spool speed from 92.9 to 93.0 
percent rated. Bleed valves closed. 


CONFIDENTIAL 




11ACA RM S54J11 


cQwwmriAL 


Outer-spool speed 
Inner- spool speed 
Turbine- Inlet temperature 
Tail- pipe temperature 
Outer-compressor-discbarge 
pressure 

Inner- compr e s sor- di s charge 
pressure 

Tail-pipe pressure 
Thrust 


Time* sec 

Figure 29. - Semilog plot of response of engine variables to change in 
exhaust -nozzle area for decrease in corrected outer-spool speed from 
65 • 9 to 63,1 percent rated and decrease in corrected inner-spool speed 
from 83.8 to 83,1 percent rated. Bleed valves open. 
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Figure 31. - Variation of outer- spool time constant with engine speed for 
change in exhaust-nozzle area. 
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Corrected outer- spool speed, percent rated 


Figure 32. - Variation of turbine- inlet- temperature rise ratio with engine 
speed for change in exhaust-nozzle area. 
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Corrected outer- spool speed, percent rated 

Figure 33. - Variation of tail- pipe-pressure rise ratio 
with engine speed for change in exhaust-nozzle area. 



Corrected outer-spool speed, N Q /-y/0^, percent rated 

Figure 34. - Variation of tail-pipe- temperature rise 
ratio with engine speed for change in exhaust -nozzle 
area. 
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Figure 35. - Variation of thrust rise ratio with engine speed 
for change in exhaust-nozzle area. 
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EXPERIMENTAL DETERMINATION OF LINEAR DYNAMICS OF 
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Abstract 


Transfer functions descriptive of the response of most engine vari- 
able? were determined from transient data that were obtained from approxi- 
mate step inputs in fuel flow and in exhaust-nozzle area. The speed 
responses of both spools to fuel flow and to turbine-inlet temperature 
appeared as identical first-order lags. Response to exhaust-nozzle area 
was characterized by a first-order lag response of the outer-spool speed, 
accompanied by virtually no change in inner-spool speed. 
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